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ABSTRACT 

New early release observations of the Cepheus E outflow and its embedded 
source, obtained with the Spitzer Space Telescope, are presented. We show the 
driving source is detected in all 4 IRAC bands, which suggests that traditional 
Class classification, although essentially correct, needs to accommodate the 
new high sensitivity infrared arrays and their ability to detected deeply em- 
bedded sources. The IRAC, MIPS 24 and 70/im new photometric points are 
consistent with a spectral energy distribution dominated by a cold, dense en- 
velope surrounding the protostar. The Cep E outflow, unlike its more famous 
cousin the HH 46/47 outflow, displays a very similar morphology in the near 
and mid-infrared wavelengths, and is detected at 24/im. The interface between 
the dense molecular gas (where Cep E lies) and less dense interstellar medium, 
is well traced by the emission at 8 and 24 /im, and is one of the most exotic 
features of the new IRAC and MIPS images. IRS observations of the North lobe 
of the flow confirm that most of the emission is due to the excitation of pure 
II2 rotational transitions arising from a relatively cold {Tex ~ 700 K) and dense 
(Nh ~ 9.6 X 10^° cm~^) molecular gas. 

Subject headings: ISM: Herbig-Haro objects — ISM: individual (Cep E) — ISM: 
jets and outflows — ISM: molecules — infrared: stars — stars: formation 

1. Introduction 

Class objects were defined a decade ago (Andre, Ward-Thompson & Barsony 1993) 
as those protostellar objects undergoing a gravitational collapse. They are characterized 
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by strong submillimeter continuum emission and relatively weak or 'non existent' infrared 
emission shortward of ~ 10/xm (Andre, Ward- Thompson, & Barsony 2000). Observationally, 
all the known Class objects show traces of active mass loss in the form of outflows (Andre, 
Ward-Thompson & Barsony 2000), and therefore, theoretical schemes to describe their evo- 
lution invoke a delicate balance between mass accretion (to form the star) and mass loss, 
to drive the outflow and remove angular momentum to sustain the star formation process 
(Smith 1998; Myers et al. 1998; Smith 2000; Proebrich et al. 2003). 

The Cepheus E (Cep E) outflow and its driving source are considered a very good 
example for studying this interface between the protostar and its mass loss, given the youth 
of its outflow, with a dynamical age of ~ 5000 yr, and the strength of its millimeter continuum 
emission, 1 Jy at 1.3mm (Chini et al. 2001; Lefloch, Eisloeffel & Lazareff 1996). At a distance 
of 730 pc, Cep E-mm is one of the brightest Class protostars known, and hkely to become 
an intermediate (~3 M©) mass star (Moro-Martin et al. 2001; Proebrich et al. 2003). 

Cep E was first identified as an outfiow based on miUimeter CO observations (Sargent 
1977; Fukui 1989), followed by near infrared and higher spatial resolution CO studies (Ho- 
dapp 1994; Eisloffel et al. 1996; Ladd & Hodapp 1997; Ladd & Howe 1997b; Noriega-Crespo, 
Garnavich & Molinari 1998). The outfiow itself is quite compact, ~1.5', and driven by the 
IRAS 23011+6126 source, also known as Cep E-mm. The outflow is deeply embedded and 
nearly invisible at visible wavelengths, with the exception of its South lobe which is breaking 
through the molecular cloud (Noriega-Crespo 1997; Devine, Rcipurth & Bally 1997; Ayala 
et al. 2000; Noriega-Crespo & Garnavich 2001) and is named Herbig-Haro (HH) object 
377. The properties of the outflow have been thoroughly analyzed in two recent studies 
(Moro-Martin et al. 2001; Smith, Proebrich & Eisloffel 2003). A comparative study of Cep 
E-mm source, in the context of other well known Class O/I sources, was recently carried by 
Proebrich et al. (2003). We will rely on these works for the interpretation of the Spitzer 
observations. 

In this Letter we present new data on the Cep E source and outfiow obtained with 
the Spitzer Space Telescope using its three instruments aboard. A brief description of the 
observations is presented in § 2, in § 3 is the analysis of the data, and in § 4 we summarize 
our results. 



2. Observations 

The Spitzer Cep E data were obtained during the in-orbit checkout and Science Verifi- 
cation period. Cep E was observed using MIPS (Rieke et al. 2004) on 2003 Nov 29, IRAC 
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(Fazio et al. 2004) on Nov 26 and IRS (Houck et al. 2004) on Nov 29. The MIPS observa- 
tions were performed in photometry mode at 24/im and fine scale 70/j,m, with 3 and 10 sec 
frame times respectively, and a total integration time of ~48 and 167 sec, accordingly. The 
IRAC observations used high dynamic range 12 sec frames with an on-source time of 108 
sec. The IRS observations are of the North section of embedded outflow (a = 23^03™13.6^, 
5=+61°42'43.5" J2000). The spectroscopic observations were tuned to obtain a high signal- 
to- noise ratio (> 50) spectrum over the 5.5 to 37//m wavelength range, based on our previous 
knowledge of the Cep E mid-IR emission (Moro- Martin et al. 2001). Ramp durations of the 
different IRS modules were adjusted to achieve this. The data from the 3 instruments were 
processed with the Spitzer Science Center Pipehne S9.5^ which provides the basic calibrated 
data, the starting point for further analysis. The IRS spectra was analyzed using ISAP ^. 



3. Results and Discussion 

3.1. The IRAC-MIPS Morphology and Photometry 

Figure 1 shows an image of Cep E in the near-IR vibrational line H2 v = 1-0 S(l) at 
2.12/im (Noriega-Crespo, Garnavich & Molinari 1998) over a 2' FOV; the position of the 
IRAS 23011+6126 (Cep E-mm) is marked with a cross. A composite image using the four 
IRAC channels is shown in Figure 2. At first sight the morphology of outfiow in the H2 near 
infrared emission is remarkably similar to that at 3.6, 4.5, 5.7 and 8//m. Some noticeable 
differences are (i) the embedded driving source is bright and undoubtly detected, (ii) there 
are some fan structures at the base of both outflow lobes, (iii) the H2 emission (green color) 
is the dominant component (see below), and (iv) the edge of the cloud , where HH 377 is 
visible, is traced by the 8//m emission (red color) and is probably related to the 7.7//m PAH 
feature. 

A similar, and a bit more ambitious because of large wavelength range coverage, is the 
comparison shown in Fig. 3, where the H2 2.12//m, IRAC 4.5//m and MIPS 24/im emission 
are shown together. The comparison is relevant because it shows the presence of the outflow 
at 24/im, despite the very bright source, and a larger extension of the warm dust at the 
edge of the cloud (red color). We performed photometric measurements of the Cep E-mm 
source using IDP3^, including PSF-fltting with apertures ~ 6" for the IRAC channels, and 
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tuning the sky background measurements to the extended emission from the outflow, which 
happens to be insignificant in relationship with the flux of the source (see below). The 
MIPS measurements at 24 and 70/xm include the outflow emission, but this contamination 
is estimated to be less than 1%. The integrated fluxes for the source are in Table 1, the 
formal uncertainties on the measurements are ~ 10% for IRAC and MIPS 24/xm. For the 
70/Lim flne scale measurement the uncertainty is likely to be higher, 20%, mostly because 
the source is near the edge of the good part of array and we are underestimating the flux 
even after the aperture correction. 

The spectral energy distribution (SED) using previous and Spitzer measurements is 
shown in Fig. 4, and is compared also with a very simple SED model. The model has been 
described in detail elsewhere (Noriega-Crespo, Garnavich & Mohnari 1998). It is based on 
an envelope made of spherically symmetry dust shells where the temperature and density are 
assumed to follow a power law and exponents of these mathematical relations are adjusted to 
flt the observations. The model in Fig. 4 has identical parameters as those in Moro-Martin 
et al. (2001), which uses silicates as the main component of the dust grains, and therefore, 
produces big absorption bands at 10 and 20yum. One of the reasons why we adopted this 
model was the fact that ISOCAM observations did not detect the source at 9.6/im (Noriega- 
Crespo, Garnavich & Molinari 1998). The model SEDs in Fig. 4 are only meant to illustrate 
what is expected from an dust envelope model and to bring into context the new IRAC 
and MIPS measurements. There is the concern that the IRAC measurements could be 
contaminated by shock emission from the outflow. In Cep E-mm source this certainly could 
be the case if a large aperture (> 6") is used (see e.g. Engelke at al. 2004, Fig. 10). Figure 5 
shows the CVF ISOCAM spectrum of Cep E-mm, the spectrum was extracted from a single 
6" X 6" pixel, and shows strong silicate, H2O and CO2 absorption features, characteristic of a 
deeply embedded source, but not traces of H2 emission from the outflow itself (Moro-Martin 
et al. 2001). We believe, therefore, that any light contribution from the outflow in the Cep 
E-mm source IRAC measurements is negligible in this case. The dust envelope SED model 
predicts a mass of ~13 Mq and a bolometric luminosity of ~30 Lq. SED models based on 
a graybody approximation (Froebrich et al. 2003) predict a factor two higher luminosity (~ 
78 Lq) and temperature (~ 35 K), although the envelope mass differs only by 25% (~8 Mq) 
with respect to our model. 

We have pointed out before (Cernicharo et al. 2000) that despite the many magnitudes of 
absorption produced by the dense envelopes surrounding Class protostars and the presence 
of deep absorption features from silicate dust and ices, that is possible to detect low mass 
Class protostars at shorter wavelengths than 10//m, given the high sensitivity of infrared 
arrays in space. This is relevant because of large surveys which are being performed by the 
Spitzer Legacy projects, aiming to provide a complete census of the low end of the mass 
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spectrum in nearby molecular clouds (Evans et al. 2003), and should be taken in to account 
when a protostar class is assigned. 

3.2. H2 Emission from Cep E North 

The IRS observations of the outflow were taken at the brightest region of the embedded 
redshifted lobe, as illustrated in Fig. 1. The spectrum of Cep E North is dominated at short 
wavelengths (<13/xm) by the pure H2 rotational hues (Fig. 6) as expected from previous 
studies (Moro-Martin et al. 2001; Smith, Froebrich & Eisloffel 2003; Noriega-Crespo 2002). 
The detection of the H2 0-0 S(0) 28.22/im and 0-0 S(l) 17.04^m lines is among the first 
ones in a low mass outflow, and demonstrates the superb sensitivity of the IRS spectrometer. 
Besides the II2 rotational lines, two atomic low excitation lines, [SI] at 25.25/im and [Sill] at 
34.8/i.m, are detected. The spectrum shows the 11.3/xm PAH band, and there are indications 
that silicate absorption at 10/xm and CO2 ice at 15.2//m are also present. 

The spectrum was dereddened using an E(B-V)= 1, and the measured H2 integrated 
fluxes were used to estimate the mean excitation temperature (assuming LTE) and the 
column density of the emitting molecular gas. We found T^^ = 686±30 K and Nh ~ 9.6 x 10^° 
cm~^, this column density corresponds to a number density of ~ 6.4 x 10^ cm~^, assuming 
an an envelope size of ~ 1000 AU (Moro-Martin et al. 2001). 

Studies of the mid and far-infrared emission from the Cep E outflow (Moro-Martin et 
al. 2001; Smith Froebrich & Eisloffel 2003) essentially agree with the scenario in which low 
velocity C-type shocks are the main excitation mechanism. Based on the relative strength 
of the H2 lines and our estimate of the gas density, we derive a mean shock velocity of 
~ 20 km s~^ for the North outflow lobe (Kaufman &: Neufeld 1996). However, the complex 
morphology of Cep E flow, which is made of multiple knots (Ladd & Hodapp 1997), may 
require as many as 20 C-type bowshocks to explain the integrated emission of their lobes 
(Smith, Froebrich & Eisloffel 2003). 

4. Summary 

We presented new observations of the Cep E outflow and its protostellar source, using 
the 3 instruments aboard the Spitzer Space Telescope. We found that: 

(1) The morphology of the outflow is remarkably similar to that of the near infrared 
observations. Considering that Cep E provides the second example of a low mass outflow 
observed with Spitzer (the other one being HH 46/47), we can speculate that outflows 
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from Class sources do not develop a "loop" or a "bubble" structure (Noriega-Crespo et al. 
2004), perhaps because they are in the initial phase of breaking through their dense placental 
envelope and futher ejection events are required to develop a cavity. 

(2) The Cep E-mm source or IRAS 23011+6126 was detected in all four IRAC chan- 
nels. Considering that the Cep E system lies at a distance of 730 pc, this suggests that 
our traditional definition of a Class source needs to accommodate the higher sensitivity 
measurements, at wavelengths shorter than 10/xm, achievable with the new generation of 
infrared arrays in space. 

(3) The IRAC and MIPS integrated fluxes of Cep E-mm source are consistent with the 
Class envelope models, although in detail our simple SEDs models do not fit the IRAC 
points. A comparison of the 5.7 and 8fim IRAC measurements with the Cep E-mm ISOCAM 
CVF spectrum shows a very good agreement. 

(4) The IRS spectroscopic observations of the redshifted North lobe show a rich spectrum 
of H2 pure rotational lines, including strong S(l) 17.04 and S(0) 28.22/xm lines. Assuming 
LTE we used these lines to estimate an excitation temperature of ~ 700K and a column 
density of Nh ~ 9.6 x 10^*^ cm~^ for a collisionally excited molecular gas. A comparison with 
C-type shock models, using the H2 line ratios and a gas density of 10^ cm~^, gives a mean 
shock velocity of 20km s~^, consistent with previous estimates. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. ANC's research was partially supported by NASA-APD Grant NRAOOOl- 
ADP-096. AAH is supported by the Spanish Program Nacional de Astronomia y Astrofisica 
under grant AYA2002-01055. 
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Table 1. Integrated mid and far- infrared fluxes" of the Cep E-mm source for the 4 IRAC 
channels, MIPS 24 and (fine scale) 70//m bands. 



Wavelength^ 3.5 


4.6 


5.7 


8 24 


70 


4.5 


29 


75 


140 3600 


30500 



^in milliJansky 

^Bandpass reference wavelength in //m 



Fig. 1. — Cep E outflow in H2 at 2.12/im . The image shows the schematic position of the 
Cep E-mm source (cross) and the Short-Low and Long-low IRS shts. The map shows a 
2' x2' field. N is up and E is to the left. 



Fig. 2. — [Color Plate] Spitzer image of the Cepheus E embedded outflow obtained using 
IRAC; the color scheme is the following: 3.6/im (blue), 4.5+5.8//m (green) and 8.0//m (red). 
The image covers a region of ~ 5' x 5', and oriented with North up and East to the left. 
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Fig. 3. — [Color Plate] Spitzer image of the Cepheus E embedded outflow obtained using 
H2 at 2.12/im (blue), IRAC Channel 2 at 4.5/im (green) and MIPS 24/im (red). The image 
covers a region of ~ 5' x 7'. North is up and East to the left 
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Fig. 4. — Spectral energy distribution of the IRAS 23011+6126 source. Data points are taken 
from Ladd & Howe 1997 and Noriega-Crespo et al. 1998 (filled diamonds), Proebrich et al. 
2003 (open triangles) and IRAC and MIPS measurements in this work (open diamonds). 
The CVF spectra of Cep E is also included (thick solid line). Two simple models are used 
for comparison, (1) a model assuming a dust opacity dominated by bare silicates, at a 
temperature of 18 K and n = 7.5 x 10^ cm"'^ (solid line), and (2) a model with n = 6.4 x 10^ 
cm~^ (broken line), 18 K and dust silicate with a thin ice mantle (Moro-Martin et. al 2001). 
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Fig. 5. — Spectrum of Cep E-mm obtained with ISOCAM Circular Variable Filter mode, 
with some strong silicate and ice absorption features (Moro-Martin et al. 2001). 
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Fig. 6. — IRS spectrum of the North lobe of the Cepheus E embedded outflow. 



